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ABSTRACT 

A numerical method of mode analysis of rapidly rotating relativistic stellar models by 
the Cowling approximation is applied to rotating neutron stars with realistic equations 
of state. For selected equations of state, eigenvalues and eigenfunctions of f-modes are 
numerically solved for stellar models from non-rotating to maximally rotating states. 

Neutral points of the lower order f-modes are determined as a function of the stellar 
rotational frequency. As in the polytropic case, we find that the bar mode can have 
neutral points for models with relatively strong gravity. The rotational frequency at the 
neutral point increases as the gravitational mass of the model becomes larger. 

As astrophysical applications of our analysis, we discuss the time scales of gravi- 
tational radiation induced instability and the possibility of the resonant excitation of 
f-modes during inspiraling motion of compact binary systems. 

Subject headings: stars:neutron — stars:oscillation — stars:rotation 
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Introduction 



In the previous paper ( Yoshida fe Eriguchi 199?| . 
hereafter YE), we determined neutral stability points 
of the Chandrasekhar- Friedman- Schutz (CFS) insta- 
bility of general relativistic rotating polytropes, by 
which non-axisymmetric oscillations of rotating stars 
are excited through the coupling with gravitational 



radiation (see Chandrasekhar 1970, Friedman fc Schutz 



1978, Friedman 1978). In the absence of viscosity 



this instability sets in at the points where eigenfre- 
quencies of the modes vanish as seen from the iner- 
tial observer at spatial infinity. Thus to determine 
the points on equilibrium sequences of rotating stars 
where the model begins to become unstable, zero fre- 
quency modes in the asymptotically inertial frame 
must be found. 

In YE the equations of state (EOS) of stellar mat- 
ter were restricted to the simple relation of polytropes 
and neutral points of the counter-rotating f-modes 
were obtained. Here the counter-rotating modes de- 
note the oscillations whose phase propagation is ret- 
rograde with respect to the stellar rotation as seen 
from the observer rotating with the star. 

In this paper we investigate oscillation modes and 
their neutral points of stability for more realistic EOS 
proposed for the neutron star matter. The investi- 
gated modes are the same ones as in YE which may 
be the most susceptible spheroidal modes to the CFS 
instability. [] 

It has been discussed that modes with the az- 
imuthal quantum number m — 3 to 5, by which 
the eigenfunction of the modes are decomposed to 
harmonics having angular (f— coordinate dependence 
~ e lmy , are the most interesting for this instability 
( Lindblom 1986| ). Also interesting is the recently dis- 
covered bar mode neutral points for rather soft EOS, 
which never appear in the Newtonian f ramework (for 
fully general rel ativistic treatments, see Btergioulas fc 
Friedman 1997 ; also see YE). Thus we will investigate 
these lower order modes in this paper. 

Concerning the neutral points to the CFS insta- 
bility, Morsink et al. (1998, hereafter MSB) investi- 
gated realistic neutron star models by applying the 
numerical method to find the exact neutral modes of 
general relativistic rotating stars developed by Ster- 



1 Recent discovery of the CFS instability of the r-mode and its 
strong effect on the stellar rotational evolution arc also interest- 
ing subject s ( A.ndersson 199s] , Andcrsson ct al. 1998, [Lindblom 
[jt al. 1998 ), mt they are beyond the scope of this paper. 



gioulas and Friedman (1997). They obtained f-mode 
neutral points for models with various masses for sev- 
eral representative EOS. Therefore we will compare 
our results with theirs. 

The counter-rotating f-modes are not only inter- 
esting in the context of the CFS instability of a single 
neutron star, but also may play an important role 
in compact binary systems because they may couple 
strongly with the tidal potential of the companion. 
The most significant effect will be the resonant exci- 
tation of the modes by the tidal force and its back 
reaction to the orbital motion of the binary system. 
We will study this subject in the last section of this 
paper. 

2. Brief Summary of the Solving Method 
2.1. Assumptions 

We assume that axisymmetric equilibrium stars are 
rotating uniformly and that the stellar matter is de- 
scribed by zero temperature EOS. Under these as- 
sumptions, equilibrium states of relativistic rotating 
stars arc obtained numerically by the KEH scheme 
( Komatsu et al. 1989 ). Linear adiabatic perturba- 
tions may be a good approximation in the present 
situation, and it is also assumed that the adiabatic 
index 7 of the perturbation coincides with the local 
adiabatic index of the equilibrium star as follows: 
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(1) 



where A means the Lagrangian perturbation of the 
corresponding variable. Eulerian perturbations of the 
metric components are totally omitted as in the pre- 
vious study (YE), i.e., the Cowling approximation is 
adopted. 

2.2. Equations of State 

There exist many candidates for EOS of real neu- 
tron stars with zero temperature. We here examine 
some of the representative EOS to cover the wide 
range of stiffness. Our choices are those of 1 ) Pand- 
haripande with hyperon (denoted by EOS B in Arnett 
fe Bowers 1977), 2) Be the- Johnson without hyperon 
flBethe fc Johnson 1974|) , 3) Bethe- Johnson (EOS C in 



2 See Nozawa et al. (1998) for recent calculation and summary 
of equilibrium models with various realistic candidates of cold 
EOS. 
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Arnett & Bowers) and 4) more recent WFF3 (Wirings 
ct $1 1988| ) joined to N V flNegel fc Vauthcrin 197 
in the low density region. In order to compare our 
result with those ol MSB, the EOS ol Pandharipande 
without hyperon (EOS A in Arnett & Bowers) is also 
employed. Extremely stiff EOS L in Arnett & Bow- 
ers is used only in Figure 8 for the comparison of the 
EOSs with a wide range of stiffness. 

2.3. Numerical Treatment 

Our numerical scheme is basically the same as that 
in YE. Perturbed quantities are assumed to behave 



as 



-27rivt+i'mtp 



where t is the killing time coor- 
dinate. A minor change that has been made is the 



introduction of a function q _ ()p/(f + p), instead of 
the Eulerian perturbation of the Emden function for 
polytropic stars. Coefficients of the perturbed equa- 
tions contain the background metric and its connec- 
tion coefficients as well as the pressure gradient and 
a function of adiabatic index like jp/(e + p). 

We have used (r x 9) — (100 x 61) grid points 
for equilibrium models where (r, 9) are the spherical 
polar coordinates. Since less number of grid point 
is used for the perturbational calculation due to the 
restrictions of the power of the computer, values for 
equilibrium states are interpolated to give values at 
the coarse grid points of our surface-fitted coordinate 
(see YE) . The interpolation is done by employing the 
cubic spline scheme in two dimensions. The results in 
this paper are obtained by using (r x 9) = (25 x 12) 
grid points in the surface-fitted coordinate. 

3. Results 

3.1. Eigenfrequencies and Eigenfunctions 

Rotational sequences of equilibrium stars can be 
obtained by fixing the central energy density e c and 
changing the rotational parameter such as the ratio 
of the polar radius to the equatorial radius in the 
meridional cross section of the star. Physical quanti- 
ties such as the gravitational mass, the T/|W| value 
and the angular momentum, etc. are calculated after 
equilibrium configurations are computed. Here T and 
W are the rotational energy and the gravitational en- 
ergy, respectively, whose definitions can be found in 
Komatsu et al. (1989) and its ratio can be considered 



to be a standard indicator of the stellar rotation. By 
changing the central density many sequences of equi- 
librium stars are obtained and a series of models with 
the same gravitational mass and a different rotational 
frequency can be chosen. For these equilibrium mod- 
els the eigenproblem is solved numerically. 

In this paper we will concentrate only on the 
counter-rotating f-modes. They are the generaliza- 
tion of the Kelvin modes of Newtonian Maclaurin 
spheroids with the counter-rotating phase velocity 
as seen from the co-rotating observer with the star 
and have indices I = m of the spheroidal harmonics. 
These particular modes may be the most suscepti- 



ble to the CFS instability (see Baumgart & Friedman 
1986) and have been mainly investigated as to this 



3 With these coefficients introduced, our system of equations is 
free from coefficients that diverge at the stellar surface for rel- 
atively stiff EOS unlike those of MSB. 



instability. 

Figures 1-7 show the dependency of eigenfre- 
qucncy v on the rotational frequency / for the spec- 
ified EOS and for the fixed gravitational mass. For 
slowly rotating models, the frequency becomes higher 
as m increases. Note that, since we are considering 
the behavior of the perturbed quantities expressed as 
^ e -2TTtvt+tmLp ^ ph ase velocity 2-Kvjm is negative 
in this case. As the rotational frequency is increased, 
the rotational dragging effect on the mode (not the 
'inertial- frame-dragging' effect in general relativity) 
works more strongly for the larger m modes and the 
order of the mode frequencies are reversed. As a re- 
sult, modes with larger m pass the neutral points ear- 
lier on the sequences. 

Near the mass-shedding limit, eigenvalues rise sharply 
with the increase of the rotational parameter. By us- 
ing more detailed data set from which these graphs 
are produced, we can have much smoother eigenfre- 
quency curves with the fixed central energy density 
and show that this behavior is also seen there. How- 
ever, the eigenfunctions of these models show sharp 
rises of the amplitudes near the surface region on the 
equatorial plane. Therefore, the rather coarse angu- 
lar resolution may prevent us from solving the nu- 
merical eigenvalue problem accurately near the mass- 
shedding limit. 

Next we show the typical behaviors of the eigen- 
frequencies of the mode of our interest. In Figure 
8 dimensionless eigenfrequency of the m = 3 mode 
is plotted against dimensionless rotational frequency 
for three of the EOSs with various stiffness, with 
the gravitational mass of the model being fixed (as 
M = IAMq). The EOS B is the softest among the 
three, and the EOS WFF3-NV, EOS L become stiffer 
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Fig. 3. — Same as Figure 1 except for M = l.OM . 



Fig. 1. — Eigenfrequencies of f- modes for the neu- 
tron stars constructed with the WFF3-NV EOS and 
M — 1.8M Q . It should be noted that both the eigen- 
frequency and the rotational frequency of the stars 
are not the angular frequencies but the ordinary fre- 
quencies. Symbols have the following meanings: 'O' 
for to = 2 mode, '+' for m = 3, '□' for m = 4 and 
'x' for to = 5. The vertical line in the right part 
is the maximum frequency which corresponds to the 
mass-shedding limit of the sequence. 




Fig. 4. — Same as Figure 1 except for the EOS B and 
M = 1.4M©. 




Fig. 2. — Same as Figure 1 except for M = 1.4M Q . 
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Fig. 5. — Same as Figure 4 except for M — 1.0M Q 
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Fig. 7. — Same as Figure 4 except for the EOS C. 




Fig. 8. — Dimensionless eigenfrequcncy of m = 3 
mode is plotted against the dimensionless rotational 
frequency. The three sequences correspond to the 
same gravitational mass (M = 1.4M©) models with 
different EOSs, WFF3-NV, B and L. 



in this order. Normalization factor of them is ^Anp, 
where p = M/V p is the averaged density with M and 
V p being the gravitational mass and the proper vol- 
ume of the equilibrium star. We can see that with 
this normalization the eigenfrequcncy is rather insen- 
sitive to the stiffness of the EOS from cases with no 
rotation to those nearly at the mass-shedding limits. 

Figure 9 displays the variation of the normalized 
eigenfrequcncy due to the gravitational mass differ- 
ence, with the EOS being fixed (as WFF3-NV). It 
is seen that increase in gravitational mass makes the 
mode frequency larger, which corresponds to the ef- 
fect of the softening of the EOS in Figure 8. This 
is reasonable since the strong self-gravity of the star 
tends to induce its density profile to be concentrated 
to the central region of it, which effectively realizes 
the configuration with the softer EOS. 

In Figures 10 and 11, we show the typical behav- 
ior of the eigenfunction q. The equilibrium models 
compared in these two figures have the same central 
density but the rotational frequencies are different. In 
both models, the function q increases monotonically 
from the center to the surface of the star along the 
radial spokes in the surface-fitted coordinate. For the 
slowly rotating model (Fig. 10), the angular depen- 
dence of the q is nearly that of the associated Leg- 
endre function, P ; m (cos#) (in this case I = m = 4), 
whereas rapid rotation tends to shift the distribution 
of the function q to the equatorial plane (Fig. 11). 

To see the 'radial' dependence of the eigenfunc- 
tion, we show the function q on the equatorial plane 
(Fig. 12). Here the same models are used as those 
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eigenfunction q : T/|W|=0.053 




Fig. 10. — An example of the eigenfunction q (see 
text) for a slowly rotating star with the WFF3- 
NV EOS. A quarter of the meridional cross section 
of the star is shown. The mode number m = 4. 
The radial coordinate distance is normalized by us- 
ing c/(47re c ) 1 / 2 where c is the velocity of light and 
e c is the central mass density of the star. The am- 
plitude of the eigenfunction is normalized such that 
the value of q at the surface point in the equatorial 
plane becomes unity. The parameters of the equi- 
librium model are: e c — 1.0 x 10 15 g/cm 3 , the ro- 
tational frequency / = 285 Hz, M = 1.19M Q and 
T/\W\ = 0.053. The eigenfrequency v = -2391 Hz. 



in Figures 10 and 11. As a 'radial' coordinate here 
we take the value x = 1 — e/e c , which roughly rep- 
resents the matter-energy distribution of the equilib- 
rium stars. As seen in this figure, the distribution 
of the function q is concentrated to the surface re- 
gion as the star rotates more rapidly. In the New- 
tonian theory, the same situation seems to improve 
the Cowling approximation for rapidly rotating stars 



(Yoshida 1997) because the role of perturbed gravi- 
tational potential becomes less important. Roughly 
speaking, less part of the stellar mass participates 
in the oscillation as the star rotates more rapidly. 
This is also likely to be the case in general relativistic 
stars, though we have no 'exact' quasi-normal modes 
to compare with. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




Fig. 11. — The eigenfunction q for a rapidly rotating 
star. The same mode as in Figure 10. The parameters 
of the equilibrium model are: e c = 1.0 x 10 15 g/cm 3 , 
/ = 1092 Hz, M = 1.43M Q and T/\W\ = 0.093. The 
eigenfrequency v = 1250 Hz. 



Fig. 12. — The eigenfunction q on the equatorial plane 
is plotted against the density coordinate x. The equi- 
librium models are the same as in Figure 10 and 11. 
Mode number m = 4. The solid curve is that for the 
model with T/|VF| = 0.053, whereas the dashed one 
is that for the model with T/\W\ = 0.093. 

In Figure 13 the differences of radial behavior in 
the equatorial plane between modes with different m 
are shown. Here the radial coordinate is that of the 
surface-fitted coordinates. It can be seen that as m 
increases, the main part of oscillation of the star shifts 
outward where the amount of the mass fraction de- 
creases. This can make the Cowling approximation 
more accurate for higher order modes (see Table 1). 

3.2. Neutral Points of the CFS Instability 

As is already remarked, we can estimate the neu- 
tral points of the CFS instability by finding zeroes of 



G 




Fig. 13. — The eigenfunction q on the equatorial plane 
with different mode numbers. The radial coordinate 
distance is that of the surface-fitted coordinate. The 
equilibrium model is the same as that in Figure 11. 
Symbols have the same meanings as in Figure 1. 

the eigenfrequencies. |^| In Table 1 we compare the 
values of T/|W| at the neutral points of the instabil- 
ity with those obtained by MSB. Here we summarize 
the tendency of the Cowling approximation in gen- 
eral relativity as follows: (1) for the same EOS and 
for the same mode number, the Cowling approxima- 
tion gives better results as the central energy density 
of the model increases; (2) for the same equilibrium 
model, it results in more accurate values for larger 
mode numbers; and (3) the Cowling approximation 
has a tendency to overestimate the stability in the 
case of relatively we ak gravity (cf. the Newtonian 
case in Yoshida 1997 ). These are qualitatively consis- 
tent with the previous results in YE. In contrast to the 
property (3), for larger central density (strong grav- 
ity) models, higher order modes seem to be underesti- 
mated in its stability by the Cowling approximation. 
As for the bar mode there seems no improvement with 
increase of the central density. 

When the EOS is fixed, we can calculate an eigen- 
frcquency of a mode for the model with a given rota- 
tional frequency / (Hz) and a given gravitational mass 
M (Mq). Then we have a neutral stability curve of 
the CFS instability for the mode in the f — M plane 
as a set of zeroes of the eigenfrequencies. Four neutral 
stability curves corresponding to four different modes 
with different EOS are shown in Figures 14 - 17. 



Bethe-Johnson (neutron) 




Fig. 14. — Neutral stability curves of the f-mode in 
the f—M plane. The gravitational mass is normalized 
by the solar mass. The EOS is that of Bethe-Johnson 
without hyperon contribution. The solid line is the 
mass-shedding limit curve. The dashed line is the ap- 
proximate maximum mass curve for a given rotational 
frequency (see text). Symbols used here are the same 
as in Figure 1. 




Fig. 15. — Same as Figure 14 except for the WFF3- 
NV EOS. 



4 Note that the neutral points here are not determined by using 
the graphs shown in the previous section. More detailed data 
sets have been used to obtain them. 



7 



Table 1: Comparison of values of T/|W| at neutral 
points with the results by Morsink et al. (1998) 




Fig. 16. — Same as Figure 14 except for the EOS 
C. The approximate maximum mass curve is omitted 
here. 




600 900 1000 1200 1400 1600 1800 2000 2200 

rotational frequency (Hz) 



Fig. 17. — Same as Figure 14 except for the EOS B. 
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Figure 14 displays those for the EOS of Bethe- 
Johnson without hyperon. The region of the left hand 
side of each curve is the stable region against the CFS 
mechanism for the corresponding mode. The solid 
line is the mass-shedding limit curve on which the 
gas on the surface in the equatorial plane rotates with 
the local Keplerian velocity. The dashed line in the 
upper region is the approximate maximum mass line 
for a given rotational frequency. The reason why we 
use the word 'approximate' is that to obtain it we do 
not constrain the rotational frequency to be constant 
but the axis ratio of the configuration to be constant 
( |Komatsu ct al. 1989| ). 

At first sight, it seems strange that for sufficiently 
rapidly rotating cases we have a model in the right 
hand side region of the mass-shedding curve. It seems 
to imply that for a given mass we have equilibrium 
models that have larger rotational frequencies than 
the mass-shedding case. This comes from the fact 
that the rotational frequency / is not a proper mea- 
sure of the rotation for extremely rapidly rotating 
stars (see appendix for its explanation). 
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4. Applications 

We here apply our results to two issues of interest 
in neutron star physics. One is the estimation of the 
time scale of the CFS instability. The other is related 
to the resonant excitation of f-modes in inspiraling 
compact binary systems which are possible targets of 
the gravitational wave detectors under construction. 

4.1. CFS Instability of Neutron Stars 

In the previous section neutral points of the CFS 
instability are determined. Then what we want to 
know next is how fast these unstable perturbations 
grow beyond these points. Unfortunately our approxi- 
mation does not provide us the answer directly. What 
is expected to happen in the real process is that the 
stellar free oscillations with non-zero frequencies cou- 
ple to gravitational wave radiation which carries its 
energy to infinity, and the frequencies inevitably have 
imaginary parts. This is the problem of the so-called 
quasi-normal mode which is well-known in black hole 
physics. To obtain complex eigenfrequencies, we must 
take the metric perturbations into account and the 
out-going wave conditions must be specified to them. 
This task is extraordinarily difficult for rapidly rotat- 
ing stars and no investigation has been accomplished 
for it yet. 

As its alternative, there have been several works 
which estimates the growth time of the instability by 
applying the gravitational radiation reaction poten- 



tial of the post-Newtonian expansion (Thorne 196E) 
to Newtonian stellar oscillations. In this approxima- 
tion, we only need to know the time dependency of 
the mass multipole of the oscillating star. 

Here we roughly estimate the growth rates of the 
instability by following Comins (1979), who investi- 
gated the secular effects of the gravitational radia- 
tion reaction and viscosity on the oscillating Maclau- 
rin spheroids. 

According to the analysis of Comins, the effect of 
gravitational radiation adds a small imaginary part to 
the oscillation frequency viewed from the co-rotating 
frame with the star £ = — 2n(v — f) as follows: 



2iG 



5S = 



A I 



4 7T/,', 1 



(m+l){m + 2)R 2 1 m+1 {2n 



(to - 1) [(2m + 1)!!] [(m - 1)/ - v] c 



2m+l 



(2) 

where c and G are the velocity of light and the grav- 
itational constant, respectively, and R\ is the equa- 



torial radius of the star. From this expression, we 
can estimate the e- folding time tq ^ for the growth of 
perturbations as follows: 
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where k(m) is defined as 
10 2m-4 ^/ 10 



101 2m+l 



fc(m) 



(m- 1) [{2m + 1)! 



4(2 7 r) 2m - 1 (m + l)(m + 2) 



(4) 

This formula is applied to the neutron star models 
of M = 1.4M with the WFF3-NV EOS (Table 2). 

It is seen that, as the rotation rate is increased, the 
lower order modes suffer stronger destabilization ef- 
fect by gravitational radiation, which originates from 
the efficiency of gravitational radiation of these lower 
order oscillations, than the higher order modes whose 
neutral points reside at the lower rotational frequen- 
cies. As a result these modes shown here have the 
same order of timescale near the mass-shedding limit. 

Table 2: Estimated timescale in units of second of 
the CFS instability for the M = 1.4M star with the 
WFF3-NV EOS. 
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TO — 3 


TO = 4 


m = 5 


818 






9 x 10 15 


952 






2 x 10 9 
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4 x 10 19 


4 x 10 7 


6 x 10 7 
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3 x 10 7 


1 x 10 6 


3 x 10 6 


1117 


1 x 10 3 


1 x 10 3 


3 x 10 3 



This estimation needs to be treated as very rough 
one because we assume that the formulae for New- 
tonian Maclaurin spheroids are applicable to the rel- 
ativistic neutron star models and because the corre- 
spondence of equilibrium quantities such as the mass 
and the equatorial radius is vague. However, qualita- 
tive behavior of the modes would be the same even if 
more refined treatment would be employed. 
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4.2. Resonant Excitation of the f-modes for 
Inspiraling Compact Binary Systems 

Inspiraling compact binary systems (neutron star- 
neutron star (NS/NS), black hole-neutron star (BH/NS) 
and black hole-black hole) are the most promising 
sources of gravitational wave for large interferometry 
gravitational wave detectors under construction such 
as LIGO and VIRGO (see e.g. ffhorne lggj and the 
references therein). These detectors will be able to ob- 
serve the inspiraling phase where components of the 
system are well approximated by the point masses. 
To extract meaningful results from the gravitational 
wave signals, it is indispensable to have sufficiently 
accurate theoretical templates of wave forms in the 
frequency range of 10 — 10 3 Hz to which these detec- 



tors are sensitive (Cutler et al. 1993). 

In this context, if at least one of the components 
is a neutron star, its internal hydrodynamical degrees 
of freedom may be a potential threat to the template 
construction. Bildsten & Cutler (1992) examined the 
problem of the 'equilibrium tide' and the issue of the 
tidal locking of the components. The tidal locking 
which causes synchronization of the spin and the or- 
bital motion seems rather unlikely to occur according 
to their result, and the theoretical template suffers 
only a negligible correction by it. Excitations of stel- 
lar oscillations of the inspiraling stars, or the 'dynam- 
ical tide', is another issue to be considered. Reiseneg- 
ger & Goldreich (1994) and Lai (1994) investigated 
resonant excitations of g-modes and their effects on 
the inspiral orbit. For slowly rotating stars g-modes 
frequencies may fall in the range as low as the or- 
bital resonant frequency. According to their results, 
however, g-modes affect the orbit negligibly, since the 
coupling between tidal potential and g-mode eigen- 
functions is rather weak. 

As we have seen, counter-rotating f-modes of neu- 
tron stars pass the neutral points viewed from the 
asymptotic incrtial frame if the star rotates suffi- 
ciently rapidly. Thus it may suffer resonant excita- 
tions during the inspiral phase|] Once the resonance 
condition is fulfilled, the f-mode is likely to be a much 
more dangerous obstacle to the construction of the 



5 As seen in Bildsten & Cutler, the tidal synchronization of in- 
spiraling components is unexpected. So the initial angular ve- 
locities of both components are preserved during the inspiral. 
We here consider the case in which at least one neutron star in 
the system has a sufficient angular velocity to suffer the orbital 
resonance. Whether such systems do exist actually or not is 
beyond our discussion here. 



theoretical template, since the f-mode eigenfunction 
couples more strongly with the tidal potential. We 
here adopt the simple oscillator model by Reiseneg- 
ger & Goldreich (1994) and examine the effect of the 
mode excitation. 

A sinusoidal external force operating on a star with 
a mass Mi whose frequency fulfills the resonance con- 
dition, v/m = n orb , where n orb is the orbital fre- 
quency, excites a mode whose energy amounts to, 



{F8tf 
8Mi 



(5) 



during the resonance time interval St, which is typi- 
cally the decay time of the orbit by gravitational radi- 
ation. The external tidal force by the companion with 
a mass Mi is estimated by the following formula: 



F 



GM1M2 (R x 



(6) 



where R\ is the stellar radius and a is the separation 
of the binary system. The factor S is the 'overlap 
integral' describing the efficiency of the tidal force on 
an eigenfunction defined by 



-gu 



drdOdip (Jj) £• P, 



(7) 



where g is the determinant of the metric of the back- 
ground spacetime, u the time component of the 4- 
velocity of the unperturbed stellar fluid, £ is the La- 
grangian displacement and P = V(r m Y™(8, ip)). 

The vibrational energy of an excited mode e is 
compared with the orbital energy decrease by gravi- 
tational radiation AE in the time interval St, 



^=- m [M 2 /M 1 ;R l] (^ rz 



(8) 



where the factor a m depends on the mode number, 
the mass ratio of the components and the stellar ra- 
dius. If this ratio is not negligible compared with 
unity, the assumption that the binary orbit evolves 
solely by gravitational radiation from the orbital mo- 
tion should be amended. If the stellar radius and 
the mass ratio of the components are assumed to be 
Ri = 10km and M2/M1 = 1, the factors a m are com- 
puted as follows: 



(9) 
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If the mass ratio is larger, say M2/M1 = 10 
(BH/NS case), they are: 



(10) 



The overlap integral S may be computed by us- 
ing the eigenfunction of the mode obtained by the 
Cowling approximation. We find that S ~ 1 for the 
bar mode and S is larger than 0.6, 0.3, and 0.2 for 
m = 3, 4, and 5 modes, respectively, for all the mod- 
els from the non-rotating state to the mass-shedding 
limit. 

4-2.1. Neutron Stars with the Prograde Rotation to 
the Orbital Motion 

When the star under consideration rotates in the 
same direction as its orbital motion, the resonance 
condition of 'counter-rotating' modes and the orbital 
motion may be fulfilled for the states beyond the neu- 
tral points of the modes. In this case, the bar mode 
may not be significant except for the stars which ro- 
tate with velocities of nearly the mass-shedding limit. 

As seen from Figs. 1-3, for higher order modes, 
the orbital frequencies at the resonances, n orb — v/m, 
may be under 10 3 Hz for the realistic EOS [], even at 
the mass-shedding limit of the star. Using the for- 
mula of e/AE above, it is observed that the higher 
modes than m = 4 may not affect the orbital evo- 
lution governed by gravitational radiation during the 
inspiral phase which will be observed by the gravi- 
tational wave antennas like LIGO/VIRGO. It is also 
seen that for systems with a large mass ratio like a 
BH/NS binary, the energy deposited in the NS vibra- 
tion is too small to affect the binary orbit. 

As for the m = 3 mode, its excitation may af- 
fect the stability and the evolution of the system if 
the star has a sufficiently large rotational frequency 
(for example, larger than 850Hz for M — I.OMq ; 
larger than 1000Hz for M = 1.4M ), since the vi- 
brational energy is no more negligible compared with 



"Here we made rough an estimation which assumes the valid- 
ity of Newtonian and slow-rotation-limit formulae used in the 
Reisenegger & Goldreich in our investigation. It is intended 
only to display that the overlap integral should be nearly order 
of unity in the low order mode cases. 

7 Here we pay attention only to the models with the WFF3-NV 
EOS. 



the gravitational radiation energy loss from the sys- 
tem. Moreover n OIb at its resonance is in the range of 
sensitivity of LIGO/VIRGO detectors, the excitation 
may be observed as the discrepancy of the signal form 
from its theoretical template. 

4-2.2. Neutron Stars with the Retrograde Rotation 
to the Orbital Motion 

If the star rotates in the retrograde direction against 
its orbital motion, the resonance condition of the 
modes and the orbit requires the rotational frequency 
of the star to be lower than that of the neutral points. 
Thus the resonance condition can drop to the fre- 
quency range of LIGO/VIRGO sensitivity windows 
even for lower m modes. 

Significant is the fact that the bar mode which cou- 
ples most strongly to the tidal potential can be res- 
onant for a wide range of rotational frequency of the 
star. For M = 1.4M Q models, stars with rotational 
frequency above 400Hz may suffer the resonant exci- 
tation on its inspiraling path in the gravitational wave 
antennas' sensitivity window (for M = l.OM case, 
this frequency may go down to 200Hz). The factor 
ai amounts to 50 for an equal mass NS/NS binary, 
and to 2 even in BH/NS cases with the mass ratio 
M2/M1 = 10. Since the overlap integral S ~ 1 for 
the bar mode, the resonant excitation would affect 
the evolution of the binary orbit significantly. 

Moreover the m = 3 mode can be excited in the 
stars rotating with frequency as low as 20Hz (for 
IAMq case). Thus theoretical templates of gravi- 
tational wave signals from compact binary systems 
containing NS with retrograde rotation should almost 
always take the resonant excitation of the counter- 
rotating f- modes into account. 
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A. Equilibrium Models with Lower Masses 
than the Mass-Shedding Limit Mass 

In Figures 14 and 15 we see a strange behavior of 
the rapidly rotating equilibrium models, i.e., the ex- 
istence of equilibrium models in the right part of the 
figure divided by the mass-shedding curve. It origi- 
nates from the fact that the rotational frequency be- 
comes an improper index for the stellar rotation for 
sufficiently rapid rotation models. 



f=1200(Hz) 

I.O(Hz) 



Fig. A18. — The constant-rotational-frequency curves 
in the central density — gravitational mass plane. 
Here the models are constructed with the EOS Bethe- 
Johnson (neutron). 




This 2-column preprint was prepared with the AAS IATjrjX 
macros v4.0. 



Central Density 

Fig. A19. — Schematic figure of constant-rotational- 
frequency curves. The solid line is the mass-shedding 
curve. The dotted and dot-short-dashed lines re- 
spectively corresponds to the constant-rotational- 
frequency curves with slow and rapid rotation. 

To see this we should know how gravitational mass 
changes with the central density of the star (Figure 
18, Figure 19). In Figure 18 gravitational mass M 
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of the equilibrium models with the same angular fre- 
quency is shown as a function of the central energy 
density. The solid line is the mass-shedding curve be- 
low which equilibrium models can be allowed to exist. 
Figure 19 is a schematic abstraction of the situation 
in Figure 18. 

While the angular velocity is small the curve with 
a constant rotational frequency / has a maximum 
in the allowed region, thus the mass of the mass- 
shedding configuration M a in Figure 19 is smaller 
than the maximum mass M x for this value of angu- 
lar frequency (M a < M < M x ). However, for the 
case with sufficiently rapidly rotation, the curve with 
constant angular frequency shows monotonically de- 
creasing behavior with no local maximum in the al- 
lowed region. In this case, the gravitational mass at 
the mass-shedding model Mt is larger than that of 
other equilibrium models with the same rotational fre- 
quency (Mb > M). This situation corresponds to the 
existence of the equilibrium model with smaller M 
than the mass-shedding limit in Figures 14 and 15. 
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